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ABSTRACT 



We present the results of a search for cool subdwarfs for which our candidates 
were drawn from a V — J reduced proper-motion diagram constructed by Salim 
& Gould (2002). Kinematic (U, V, and W) and self-consistent stellar parameters 
(T cff , log g, [Fe/H], and £ t ) are derived for 91 candidate subdwarfs based on high 
resolution spectra. The observed stars span 3900K < T eff < 6200K and —2.63 < 
[Fe/H] < 0.25 including only 3 giants (log g < 4.0). Of the sample, 77 stars 
have MgH lines present in their spectra. With more than 56% of our candidate 
subdwarfs having [Fe/H]< —1.5, we show that the V — J reduced proper-motion 
diagram readily identifies metal-poor stars. 
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1. Introduction 

The chemical history of our Galaxy has been the target of numerous theoretical and 
observational studies for many decades. Much progress has taken place where the large and 
increasing body of observed stellar abundances (e.g., Edvardsson et al. 1993; McWilliam 
et al. 1995; Ryan et al. 1996; Reddy et al. 2003) is complemented by the ongoing refinement 
of Galactic chemical evolution models (e.g. Timmes et al. 1995; Goswami & Prantzos 2000; 
Alibes et al. 2001). The ultimate objective of Galactic chemical evolution studies is to be 
able to predict accurately elemental and isotopic abundances for any given location and 
time. The spectra of cool stars show molecular bands from which measurements of isotopic 
abundance ratios can be made (e.g., Mg from MgH, Ti from TiO). These isotopic ratios offer 
a unique opportunity to test models of Galactic chemical evolution since different stellar sites 
and nucleosynthetic processes may be responsible for production of the individual isotopes. 
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Comparing the observed and predicted evolution of isotopic abundances will test aspects of 
the contributions of stellar nucleosynthesis to Galactic chemical evolution in a way that is just 
not possible with elemental abundances. In order to measure isotopic ratios from molecular 
lines, cool subdwarfs (unevolved metal-poor stars that are located below the solar-metallicity 
main sequence in color-magnitude diagrams) must first be found. 

Previous studies such as Ryan (1989) and Carney et al. (1994) drew candidate subdwarfs 
from proper-motion catalogs and were successful in identifying large numbers of metal-poor 
stars. Subsequent detailed abundance analyses were carried out on these targets. Both 
studies focused upon warmer targets, neglecting cooler subdwarfs. Objective prism surveys 
which identify metal-poor stars by the weakening of Ca II H and K have also been successful 
in identifying considerable numbers of metal-poor stars (e.g., Bond 1970; Beers et al. 1992; 
Christlieb 2000). Without a reliable temperature estimate, a cool metal-poor star will have 
Ca II lines comparable to a warmer metal-rich star thereby biasing the sample towards 
warmer stars. Numerous stars of spectral type M and later have recently been identified (e.g., 
Delfosse et al. 1999, Hawley et al. 2002, Reid & Cruz 2002). However, the overwhelming 
majority of these cool stars are solar metallicity. Only a handful of metal-poor stars with 
spectral type M or later have been found (e.g., Gizis 1997, Schweitzer et al. 1999, Lepine 
et al. 2003, and Burgasser et al. 2003). Our target stars are cooler than those found by Ryan 
(1989) and warmer than those found by Gizis (1997). 

In Paper I (Yong & Lambert 2003), we presented the results from our first search for cool 
subdwarfs. From these subdwarfs we have measured the evolution of the Mg isotopic ratios, 
extending the Gay & Lambert (2000) measurements to [Fe/H]=— 2.5 (Yong & Lambert in 
prep). In this paper, we present the results from our second search for cool subdwarfs, 
self-consistent stellar (T eff , log g, [Fe/H], and £ t ) and kinematic (U, V, W) parameters for a 
further 91 candidates. 



2. Selection criteria 

In Paper I, our goal was to identify subdwarfs with T e fr < 4700K and [Fe/H] < —1.5, 
that is, metal-poor stars with MgH molecular features. Following Ryan (1989), we drew 
candidates from the New Luyten Two- Tenths (NLTT) catalog of stars with annual proper- 
notions in excess of O'.'18/year (Luyten 1979,1980; Luyten & Hughes 1980). We used the 
reduced proper-motion criterion to ensure that the selected stars had large transverse veloc- 
ities relative to the local standard of rest. Since halo stars are on plunging orbits, they are 
over-represented in proper-motion catalogs. Hertzsprung first introduced the concept of the 
reduced proper- mot ion, defined as H = V + 5 log \i + 5, where V is the apparent magnitude 
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and \x is the proper-motion in arcsec/year. The reduced proper-motion can also be expressed 
as H = My + 51ogz/y — 3.37 where My is the absolute magnitude and v T is the transverse 
velocity in km/s. At a given spectral type, a reduced proper-motion constraint rejects stars 
whose transverse velocities fall below a chosen value. Subdwarfs are less luminous than disk 
dwarfs at a given color and have larger transverse velocities. Both of these factors lead to an 
increase in the reduced proper-motion of the subdwarf population relative to the disk dwarf 
population. While our original search, Paper I, was successful in finding these subdwarfs, 
our efficiency was less than ideal. Of our 134 candidates, only 11 stars had [Fe/H] < —1.5 
and MgH lines. 

For this search, we again utilized the reduced proper-motion criterion to identify metal- 
poor stars drawn from the NLTT catalog. Instead of relying upon Luyten's original photom- 
etry, we made use of the excellent work by Salim & Gould (2002, 2003) and Gould & Salim 
(2003). Salim & Gould cross-referenced the NLTT catalog with 2MASS (Skrutskie et al. 
1997) and USNO-A (Monet et al. 1998) to obtain optical and infrared photometry. Proper- 
motions were determined from the difference in positions between the two catalogs. Salim 
& Gould (2002) then constructed an optical-infrared reduced proper-motion diagram from 
which they identify "distinct tracks for white dwarfs, subdwarfs, and main-sequence stars". 
A reduced proper-motion diagram based on the original Luyten photometry does not show 
the expected separation between subdwarfs and main-sequence stars. Salim & Gould suggest 
that the reduced proper-motion diagram based on the original Luyten photometry suffers 
from a "short color baseline and large errors". In Paper I, we confirmed the Ryan (1989) 
statement that the photometric "values tabulated in the NLTT catalog must be regarded as 
approximate only". Another advantage of the optical-infrared reduced proper-motion dia- 
gram is that the colors provide an accurate estimate of the effective temperature. In Figure 
1, we show a reduced proper- motion diagram identical to that of Salim & Gould with one ex- 
ception. Rather than using V + 5 log /z, we use J + 5 log /x as the reduced proper- motion. The 
rationale is that the 2MASS photometry is more accurate than USNO. Using photographic 
plates, the principal objective of the USNO catalog is astrometry while the photometry is a 
by-product. USNO-A provides B v and Ru magnitudes which are converted to V according 
to the relation given in Salim & Gould (2000) V = R v + 0.23 + 0.32(5 - R) v . Salim & 
Gould (2003) show that the uncertainties in V are around 0.25 mag while uncertainties in 
J are around 0.02 mag. However, there are no discernable differences between the original 
(V + 5 log fi) reduced proper-motion diagram and ours ( J + 5 log //). The locations of white 
dwarfs, subdwarfs, and main-sequence stars (solar metallicity) can be easily estimated. 

Of the stars we observed in Paper I, a subset had optical and infrared photometry 
from Salim & Gould (2002, 2003). After constructing a reduced proper-motion diagram, 
we overplotted the stars we observed in Paper I. It was clear that the metal-poor stars 
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were separated from the metal-rich stars. For this search, we selected stars with Hj > 
4 + 2 x (V - J) where Hj = J + 51og/i with V and J taken from Salim & Gould (2002, 
2003). This cutoff was the approximate line dividing the subdwarfs from the main sequence 
stars studied in Paper I. To ensure we observed cool stars, we generally observed targets 
which had 1.4 < V — J < 2.6. The magnitude limit was V < 13.9, the limit of the 
McDonald 2.7m telescope. We also imposed a declination limit —40° < 5 < 90° appropriate 
for McDonald observatory. This resulted in a list of 470 candidate subdwarfs. 

In Paper I, the most metal-poor star we observed which showed MgH lines was G39-36, 
T e ff=4200K and [Fe/H]=— 2.5. Not only did this star have a large reduced proper-motion, 
it was also marked by a large ultraviolet excess, S(U — B) 06 = 0.25 (Sandage & Kowal 
1986). Roman (1955) found a correlation between the ultraviolet excess and the strength 
of metallic lines in high velocity stars. Schwarzschild et al. (1955) showed that a lower 
metal abundance and the corresponding reduction in line blanketing was responsible for the 
ultraviolet excess (see also Carney 1979 and references within). Therefore, we also observed 
cool stars with large ultraviolet excesses taken from Sandage & Kowal (1986). Finally, 
bright solar metallicity stars expected to have strong MgH lines were also added to the list 
of candidates to be observed. These bright stars were taken from the Hipparcos catalog 
(ESA 1997) and from Carney et al. (1994). In total we observed 91 targets; 58 selected from 
the V — J reduced proper-motion diagram, 24 bright targets, and 9 from (Sandage & Kowal 
1986). 



3. Observations and data reduction 

Table 1 contains the list of candidates observed at McDonald Observatory on the 2.7m 
Harlan J. Smith telescope between September 2002 and February 2003. The data were 
obtained using the cross-dispersed echelle spectrometer (Tull et al. 1995) at the coude f/32.5 
focus with a resolving power of either 30,000 or 60,000. The detector was a Tektronix 
CCD with 24 fxm 2 pixels in a 2048 x 2048 format. Using this setting, the spectral coverage 
was from 3800A to 8900A with gaps between the orders beyond 5800A. These observations 
incorporated the MgH A — X lines near 5140A. Visual inspection of the data indicated the 
presence or absence of MgH molecular features. Of our sample of 91 stars, 14 stars did not 
show MgH lines and 77 stars showed MgH lines. Numerous Fe I and Fe II lines were available 
in the observed spectra for spectroscopic determination of the stellar parameters. For each 
star, exposure times were generally 20-30 minutes and multiple exposures were taken and 
co-added to increase the signal-to-noise ratio when necessary. Although varying from star 
to star, the typical signal-to-noise ratio of the extracted one dimensional spectra was 60 per 
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pixel at 6500A. One dimensional wavelength calibrated normalized spectra were extracted in 
the standard way using the IRAF 1 package of programs. Equivalent widths were measured 
using IRAF where in general Gaussian profiles were fit to the observed profile. 

4. Analysis 

4.1. Deriving stellar parameters 

We determined the stellar parameters using the same method described in Paper I. The 
local thermodynamic equilibrium (LTE) stellar line analysis program Moog (Sneden 1973) 
was used in combination with the adopted model atmosphere. For log g > 3.5, we used the 
Nextgen model atmosphere grid for low mass stars computed by Hauschildt et al. (1999). 
For giants with log g < 3.5 we used LTE model atmospheres computed by Kurucz (1993). 
In both cases we interpolated within the grid when necessary. Equivalent widths (EWs) of 
35 Fe I lines and 5 Fe II lines were measured. The gf values of the lines were taken from 
Lambert et al. (1996) and from a compilation by R.E. Luck (1993, private communication). 
The effective temperature (T e g) was set by insisting that the abundances of individual Fe 
lines be independent of lower excitation potential. Ideally, the Alonso et al. (1996, 1999) 
T e g:[Fe/H]: color relations based on the infrared flux method would be used to obtain values 
for T e ff . In particular, the (V — K) color provides the greatest accuracy. Unfortunately, for 
[Fe/H]< -1.5, the Alonso et al. (1996) (V - K) relation is only applicable above T cff ~4600K 
and Salim & Gould (2003) have shown that errors in V are 0.25 mag. The microturbulence 
(£ t ) was set by requiring that the abundances of individual Fe lines show no trend against 
equivalent width. By forcing agreement between the Fe abundance derived from neutral and 
ionized lines, the gravity (log g) was fixed. This process required iteration until a consistent 
set of parameters were obtained (see Figure 2) from which the Fe abundance was determined 
from the mean of all Fe lines. 

Only weak lines, EWs < 90 mA, were used in the analysis if possible. We observed the 
solar spectrum at R=60,000 to check our analysis techniques. We measured 30 Fe I lines and 
7 Fe II lines and compared our equivalent widths with the Grevesse & Sauval (1999) values. 
Our equivalent widths were larger by a mean value of 3.7 mA with a standard deviation 
of 2.6 mA. We note that our equivalent widths were measured from a disk integrated solar 
spectrum whereas the Grevesse & Sauval (1999) equivalent widths were measured at disk 



1 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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center from the Delbouille et al. (1973) solar atlas. Using our equivalent widths and a 
Nextgen model atmosphere, we derived a solar abundance of loge(Fe) = 7.54. Considering 
the Grevesse & Sauval (1999) value of loge(Fe)= 7.50 ± 0.05 derived from their empirical 
model solar atmosphere, we adopted loge(Fe)= 7.52 as the solar Fe abundance for this study. 

The derived model parameters have uncertainties which were estimated in the following 
way. We varied T c g until the trends between lower excitation potential and abundance 
were unacceptable. Likewise, we changed the microturbulence until the trends between 
equivalent width and abundance were poor. The gravity was adjusted until the difference 
between the abundance from Fe I and Fe II lines was greater than the standard deviation of 
the Fe abundance derived from Fe I lines (typically 0.1—0.15 dex). Estimated uncertainties 
in the model parameters are <5T eff =150K, Slog g=0.3 dex, 5£ t =0.3 km s -1 , and 5[Fe/H]=0.2 
dex. Though the targets were fainter than those in Paper I, we increased the number of 
exposures to ensure the signal-to-noise ratio was similar. And so the uncertainties in the 
model parameters are identical to those in Paper I. 

4.2. Comparison with literature 

A search on SIMBAD indicated that 10 stars had previously determined temperatures 
and 17 stars had previously determined metallicities. In Table 2, we compare our values 
with those found in the literature. Combining our data with Paper I, we find a mean 
offset (T c fj(Yong & Lambert)— T e g (literature)) = —16 K with a standard deviation of 115K 
(see Figure 3). Again combining our data with Paper I, the mean offset is ([Fe/H] (Yong 
& Lambert) — [Fe/H] (literature)) = —0.06 dex with a standard deviation of 0.38 dex (see 
Figure 4). There is a reasonable agreement between the stellar parameters derived in this 
study (and Paper I) and the values found in a variety of sources in the literature. 

4.3. Self-consistency check 

Molecular lines will not be visible in the spectra of stars with sufficiently high temper- 
atures and/or low abundances. A simple self-consistency test of our stellar parameters is to 
predict the limit in temperature and metallicity at which molecular lines will not be observed. 
Once a detection limit is obtained, we can then check the stellar parameters to determine 
whether molecular lines should or should not be present. We synthesized representative MgH 
and TiO lines assuming a resolving power of 60,000 and log g = 4.5. The MgH molecular 
data were taken from Gay & Lambert (2000) while the TiO molecular data were taken from 
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J0rgensen (1994). Values of [Ti/Fe] and [Mg/Fe] in accord with observed trends summarized 
in Alibes et al. (2001) were assumed, i.e., [Mg/Fe]=0.4 at [Fe/H]=— 1.5. At a given value 
for T c g, we decreased the metallicity until the molecular features reached a depth of ~5% 
relative to the continuum. We considered this value of T c g- and [Fe/H] to be our detection 
limit. In Figure 5, we plot our stars in temperature-metallicity space imposing our limits of 
detection for MgH and TiO. At a fixed metallicity, MgH is detectable to considerably higher 
values of T cfr than TiO. Despite TiO (6.9 eV) having a higher dissociation energy than MgH 
(1.34 eV), the overwhelming H abundance ensures the continuing presence of MgH to higher 
temperatures. 

All 14 stars that show neither MgH nor TiO lines occupy a region in temperature- 
metallicity space where, for the given stellar parameters, we would not expect to see molec- 
ular lines. All 77 stars that show MgH lines lie in the region where, for the given stellar 
parameters, we do expect to see molecular lines. Our approximate limits of detection as- 
sume a resolving power of 60,000 and log g = 4.5. For objects observed at a lower resolving 
power, we expect the detection limit would move to lower values of T c ff at a fixed metallicity. 
Likewise, the detection limit for giants would be shifted to lower values of T c g at a fixed 
metallicity. The vast majority of the stars observed in Paper I also lie in the expected re- 
gions. That is, stars with MgH lines lie to the left of the detection limit while stars without 
MgH lines fall to the right of the detection limit in Figure 5. We note that the handful 
of stars that do not conform to the predictions are observed at lower resolving powers or 
are giants. Through the synthesis of representative molecular lines, we have performed a 
self-consistency check validating our derived stellar parameters. 



5. Discussion 
5.1. Kinematics 

We used the reduced proper-motion constraint to select stars kinematically distinct from 
the thin disk. This selection criterion was used in order to observe stars on halo-like orbits 
with halo-like metallicities. The first check we can implement is to determine whether we 
successfully targeted stars with kinematics unlike thin disk stars. Following the Johnson & 
Soderblom (1987) recipe, we calculated the Galactic space- velocity components U (positive 
towards the Galactic center), V (positive in the direction of Galactic rotation), and W 
(positive towards the north Galactic pole) along with the associated uncertainties au, <?v, 
and a w (see Table 1). For the solar motion with respect to the local standard of rest (LSR), 
we assumed the Dehnen & Binney (1998) values (+10, +5, +7) km s" 1 in (U,V,W). Note that 
in Paper I, we incorrectly applied a correction of —10 km s" 1 in U. Therefore, +20 km s -1 



- 8- 



must be added to our values of U in Paper I to obtain U lsr- In the absence of Hipparcos 
parallaxes, spectroscopic parallaxes were determined by using the derived model parameters 
and the Y 2 isochrones (Yi et al. 2001). In Figure 6 we plot Ulsr, Vlsr, an d ^lsr versus 
[Fe/H] including the stars from Paper I. As in Paper I, we identify stars that lag the LSR, 
V < —50 km s— 1, as likely members of the thick disk or halo. As expected, our reduced 
proper-motion constraint enabled us to select stars belonging to populations kinematically 
distinct from the thin disk. 



5.2. Metallicity 

Our goal was to find stars with [Fe/H]< —1.5. The second check we can conduct is to 
verify the success of our selection criteria in targeting metal-poor stars. In Figure 7 we plot 
the number of stars versus metallicity for this study, Paper I, Carney et al. (1994), and the 
Ryan & Norris (1991) sample. In Paper I, we noted that our distribution was similar to the 
Carney sample. We also mentioned that the Ryan sample had a peak at lower metallicity 
since these stars were observations of metal-poor candidates identified by Ryan (1989) as 
having an ultraviolet excess corresponding to [Fe/H]< —1.2, that is, 5(U — B) . 6 > 0.2. In 
this study, we were far more successful in identifying metal-poor stars with 38 out of 91 stars 
being more metal-poor than [Fe/H]= —1.5. As a comparison, in Paper I we had 27 of 134 
stars with [Fe/H]< —1.5. In this study, 32 stars were more metal-poor than [Fe/H]= —1.5 
with MgH lines present in their spectra. As a comparison, 11 stars in Paper I were more 
metal-poor than [Fe/H]= —1.5 with MgH lines present in their spectra. Clearly the selection 
technique utilized in this paper was more successful than the technique applied in Paper I. 
That is, the Salim & Gould (2002) optical-infrared reduced proper-motion diagram separates 
subdwarfs from main-sequence stars more effectively than a reduced proper-motion diagram 
based on the original Luyten photometry. 

Let us consider only those stars with 2MASS and USNO photometry, i.e., only the stars 
in Figure 8. After combining this study with the stars in Paper I, we find that 56% (43 of 
76) of stars that lie below the cutoff have [Fe/H]< —1.5. While this is far from the ultra 
metal-poor stars currently being discovered by other programs, [Fe/H]< —1.5 is a metallicity 
regime in which the Mg isotope ratios have not been explored. In this metallicity range, only 
cool subdwarfs will show MgH lines. We also note that beyond J + 5logfi = 10, all of the 
observed stars are more metal-poor than [Fe/H]= —1.5. 
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5.3. Temperature 

Our goal was to find stars with T c g<4700K. More specifically, ideal targets were stars 
with [Fe/H]< —1.5 with MgH molecular lines. Of the 91 stars presented in this paper, 
38 had T eff <4700K and given the uncertainties in our derived effective temperatures, we 
note that 48 of the 91 stars had T cff <4800K. As a comparison, in Paper I 44 of the 134 
stars had T cfr <4700K and 69 of the 134 stars had T cff <4800K. In Figure 9, we plot the 
number of stars versus T cff for this study, Paper I, and the Carney et al. (1994) sample. Our 
temperature distributions are similar between this study and Paper I. Both of our searches 
have focused upon stars considerably cooler than the Carney study which highlights the 
different temperature regimes of interest. Importantly, 77 of the 91 stars presented here 
have MgH lines. In Paper I, 100 of the 134 stars had MgH lines. 

In Paper I, we confirmed Ryan's findings that the Luyten colors were inaccurate. Firstly, 
the magnitudes were published to only 0.1 mag. Secondly, we chose stars assigned by Luyten 
to color classes g-k, k, or k-m. These stars spanned 4200K < T c g< 6400K yet the Luyten 
colors for this range of stars were m pg — mji=0.9, 1.1, or 1.3. We now comment on the 
colors in the Salim & Gould (2002) study. In Figure 10, we plot the V — J color against 
our derived temperatures (this study and Paper I) separated into four metallicity bins. This 
Figure shows that a fair estimate of the effective temperature can be obtained from the 

V — J color. Alonso et al. (1996, 1999) derived T eS :[Fe/Hj: color relations for dwarfs and 
giants using the infrared flux method. They have shown that the V — K color index provides 
the most accurate temperature estimates. (We refrain from using the infrared flux method 
calibrations due to the uncertainties in V .) In Figure 11, we plot the V — K color against our 
derived temperatures and find that the dispersion about the mean is slightly less than for 

V — J. The K magnitudes were also from 2MASS and presented by Salim & Gould (2002). 
This reduced scatter in the V — K versus T c g- plot suggests that a reduced proper-motion 
diagram featuring V — K versus K + 5 log \i may be a power tool when the metal-poor targets 
need to span a small temperature range. In this Figure we also plot the Alonso et al. (1996) 
(V — K):T eS relation for dwarfs. We find that this relation reasonably matches our data 
though at the metal-poor end, our data extend to cooler temperatures than the Alonso et al. 
(1996) relations. In both the V — J and V — K plots (Figures 10 and 11), we suspect that 
the dispersion is primarily due to errors in V. Earlier we mentioned that the errors in V 
were around 0.25 mag while the errors in J (and presumably K) were only 0.02 mag (Salim 
& Gould 2003). While reddening may also be a factor, we note that these cool dwarfs lie 
within lOOpc of the sun and errors in V — J and V — K would overwhelm any reddening 
effects. 

Based on spectrum synthesis, Cottrell (1978) showed that MgH lines in cool stars can 
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be a useful metal abundance discriminant. In cool stars, the strength of the MgH lines 
does not strongly depend on metallicity, unlike atomic lines. For sufficiently cool stars, a 
decrease in the metal abundance will weaken the atomic lines whilst the MgH lines remain 
strong. In Figure 12, we show the spectra of 4 stars highlighting some MgH lines and a 
strong atomic line. In 2 of the spectra the "normal" situation is shown in which the MgH 
lines are dwarfed by the strong atomic line. In the other 2 spectra, Cottrell's prediction 
is verified since the MgH lines are comparable in strength to the atomic line. That is, the 
atomic line has weakened while the MgH lines have maintained their strength. Note that all 
stars shown in Figure 12 are very cool. If these stars were solar metallicity, Figure 5 shows 
that the spectra of these stars would be swamped with TiO lines. 



6. Concluding remarks 

We present stellar parameters for 91 candidate subdwarfs selected primarily by their 
large reduced proper-motions. Our goal was to identify metal-poor stars [Fe/H]< —1.5 with 
MgH lines. This is the 2nd search that we have conducted for these cool subdwarfs. In our 
first search, we drew candidates from a reduced proper-motion diagram based on Luyten's 
original photometry. While we were able to find cool subdwarfs, our efficiency was low (11 
of 134 stars had [Fe/H]< —1.5 and MgH lines). Errors in the Luyten photometry resulted in 
a reduced proper-motion diagram without cleanly separated populations of subdwarfs and 
solar-metallicity main-sequence stars. In this search, we made use of the 2MASS infrared and 
USNO optical photometry presented by Salim & Gould (2002). The reduced proper- motion 
diagram based on this photometry showed distinct populations of subdwarfs and main- 
sequence stars. Spectroscopic observations of these candidate subdwarfs showed that the 
majority were metal-poor stars. Of the candidate subdwarfs selected from their location in 
the optical-infrared reduced proper-motion diagram, we have shown that 56% have [Fe/H]< 
— 1.5. In Paper I, only 20% of candidates had [Fe/H]< —1.5. Further, we verified that the 
V — J color index was well correlated with T e g such that cool stars suspected of having MgH 
lines could be targeted. We intend to make further observations to identify cool subdwarfs. 
In this search, we found a further 32 stars with [Fe/H]< —1.5 and MgH lines. We intend to 
re-observe candidates from this study with higher resolving power and higher signal-to-noise 
ratio in order to measure the Mg isotope ratios to study Galactic chemical evolution. 

We are grateful to Samir Salim and Andy Gould for providing a list of candidate sub- 
dwarfs prior to publication. We thank G. Fritz Benedict for valuable discussions and the 
referee, Andy Gould, for helpful comments. We acknowledge support from the Robert A. 
Welch Foundation of Houston, Texas. This research has made use of the SIMBAD database, 
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Fig. 1. — V — J reduced proper-motion diagram. The points are NLTT stars with 2MASS 
and USNO-A identifications make by Salim & Gould (2002, 2003). Main-sequence stars 
(MS) are separated from the subdwarfs (SD) by the solid line. White dwarfs (WD) are 
separated from SDs and MS stars by the dashed line. 
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Fig. 2. — Determination of stellar parameters T e ff, log g, and £ t using excitation and ioniza- 
tion equilibria for G266-8. In the top panel, the lower excitation potential (LEP)-abundance 
relation is used to set T e s- In the middle panel, the reduced equivalent width (W\/X)- 
abundance relation is used to determine In the bottom panel, the abundances of Fe I 
(asterisks) and Fe II (filled circles) are used to fix log g. In all panels the line represents the 
linear least squares fit to the data. 
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Fig. 3. — T cff comparisons between Yong & Lambert (this study and Paper I) and various 
literature studies. The open circles represent the data from Paper I while the filled circles 
are stars from this study. The solid line represents the line of equality. 
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Fig. 4. — [Fe/H] comparisons between Yong & Lambert (this study and Paper I) and various 
literature studies. The open circles represent the data from Paper I while the filled circles 
are stars from this study. The solid line represents the line of equality. 



-18- 




3500 4000 4500 5000 5500 6000 6500 



Fig. 5. — Limit of detection for MgH and TiO features based on spectral synthesis of rep- 
resentative lines assuming log g = 4.5 and R=60,000. The filled circles are from this study 
while the open circles are from Paper I. With the exception of a few stars from Paper I, all 
stars with MgH lines lie to the left of the MgH detection limit. Also, all stars without MgH 
lines lie to the right of the MgH detection limit. The few stars that do not lie where they 
ought to are giants or were observed at lower resolving powers. 
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Fig. 6. — Galactic space- velocity U, V, and W (relative to the LSR) vs [Fe/H]. The open 
circles represent the data from Paper I while the filled circles are stars from this study. 
Only stars which have ^j^f^*j^| < 0.7 are shown. A considerable fraction of the sample 
noticeably lag the LSR (V < —50 km/s). 
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Fig. 7. — Number of stars versus metallicity for Paper I (top panel), this study (2nd panel), 
Paper I and this study (3rd panel), only the stars falling below the cutoff defined in Figure 1 
(4th panel), the Carney et al. (1994) study (5th panel), and the Ryan & Norris (1991) study 
(bottom panel). The 4th panel shows most clearly that our selection technique efficiently 
selects metal-poor stars. 
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Fig. 8. — V — J reduced proper motion diagram based on the Salim & Gould (2002) pho- 
tometry. In each panel, the small dots are stars in the NLTT catalog with 2MASS and 
USNO photometry. We overplot (large circles) the locations of stars we previously observed 
for which we were able to derive metallicities either in this study or in Paper I. In the upper 
left panel the large circles are stars with [Fe/H]< —1.5, in the upper right panel are stars 
with —1.5 <[Fe/H]< —1.0, in the lower left panel are stars with —1.0 <[Fe/H]< —0.5, and 
in the lower right panel are stars with [Fe/H]> —0.5. The dotted and dashed lines (based 
on Samir & Gould) define approximate boundaries between main sequence stars (MS), sub- 
dwarfs (SD), and white dwarfs (WD). In the upper left panel, the stars with [Fe/H]< —2.0 
(open circles) overlap with the [Fe/H]< —1.5 stars. 
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Fig. 9. — Number of stars versus effective temperature for Paper I (top panel), this study 
(2nd panel), Paper I and this study (3rd panel), and the Carney et al. (1994) study (bottom 
panel). The filled histograms represent the distribution of stars with [Fe/H]< —1.5. Clearly 
our samples peak at lower values of T e g- than the Carney sample. The 2nd panel shows that 
our selection technique is very efficient at identifying cool stars with [Fe/H]< —1.5. 
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Fig. 10. — V — J color from Salim & Gould (2002) versus effective temperature derived 
from this study and Paper I. In the upper left panel are stars with [Fe/H]< —1.5, in the 
upper right panel are stars with —1.5 <[Fe/H]< —1.0, in the lower left panel are stars with 
— 1.0 <[Fe/H]< —0.5, and in the lower right panel are stars with [Fe/H]> —0.5. Given a 
value of V — J, the corresponding T c fj can be determined with reasonable accuracy. 
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Fig. 11. — V — K color from Salim & Gould (2002) versus effective temperature derived 
from this study and Paper I. In the upper left panel are stars with [Fe/H]< —1.5, in the 
upper right panel are stars with —1.5 <[Fe/H]< —1.0, in the lower left panel are stars with 
— 1.0 <[Fe/H]< —0.5, and in the lower right panel are stars with [Fe/H]> —0.5. In each 
panel, the solid line is the (V — K):T e g relation from Alonso et al. (1996) for dwarfs. The 
dotted line is the extrapolation of the (V — K):T ef[ relation. 
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Fig. 12. — Observed spectra of 4 different cool subdwarfs. The effective temperature and 
metallicity of each star is shown. In each panel, the 3 arrows highlight representative MgH 
lines while the asterisk highlights an atomic line (Fe and Ni). In the upper panels, the 3 
MgH lines are dwarfed by the strong atomic line. In the lower panels, the 3 MgH lines are 
of comparable strength to the atomic line. 
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Table 1. Basic data and derived parameters for objects 
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-1.37 


4.5 


0.4 


-94 


33 


-345 


94 


-61 


5 


Y 


Samir02 


LP 319-42 


115232 


273051 


4600 


-1.8 


4.5 


0.4 


-184 


45 


-203 


26 


7 


50 


Y 


Samir02 


LHS 343 


125624 


154143 


3900 


-2.63 


4 


0.5 


190 


60 


-625 


230 


59 


44 


Y 


Samir02 



-27- 



Table 1 — Continued 



Name R.a a Dcc b T cff [Fc/H] log g £ t Ulsr &u Vlsr. °"V Wlsr a W Molecular Source 





2000 


2000 


(K) 




(cm s 2 ) 








(km/s) 








Features 




ROSS 471 


132041 


-122606 


4600 


-2.13 


4.5 


0.5 


-88 


37 


-97 


43 


-15 


6 


Y 


Samir02 


G 63-40 


133648 


190513 


4600 


-1.95 


4.5 


0.3 


-160 


55 


-351 


134 


-58 


12 


Y 


Samir02 


NLTT 35929 


135935 


224757 


5200 


-0.57 


4.5 


0.5 


-38 


21 


-137 


77 


-13 


14 


Y 


Samir02 


G 136-67 


150816 


150948 


5000 


-0.5 


4.5 


0.6 


-88 


22 


-65 


18 


43 


18 


Y 


Samir02 


G 181-19 


165910 


345200 


4800 


-1.66 


4.5 


0.3 


-18 


13 


-324 


71 


129 


70 


Y 


Samir02 


G 140-16 


175248 


143837 


6000 


-0.32 


3.75 


1.4 


842 


754 


-835 


689 


-471 


390 


N 


Samir02 


LP 393-3 


193307 


235951 


4300 


-1.92 


4.5 


0.7 


_4 


21 


-260 


12 


-67 


12 


Y 


Samir02 


G 23-7 


193544 


14342 


5400 


-0.68 


4.5 


0.6 


33 


20 


-170 


26 


-24 


8 


Y 


Sandage86 


LP 813-13 


194343 


-155524 


6000 


-1.42 


4.5 


1 


138 


31 


-292 


83 


-75 


22 


N 


Samir02 


HIP 97174 


194503 


291943 


4700 


-0.83 


4.5 


0.3 


-105 


6 


-57 


3 


12 


2 


Y 


Hipp 


NLTT 48700 


200516 


311210 


4800 


-1.9 


4.5 


0.2 


-275 


4 


-232 


50 


40 


37 


Y 


Samir02 


G 231-39 


211833 


522336 


4600 


-1.3 


4.5 


0.3 


260 


58 


-226 


4 


9 


6 


Y 


Samir02 


NLTT 51106 


212214 


-270217 


4800 


-1.94 


4.5 


0.4 


-223 


90 


50 


18 


-182 


67 


Y 


Samir02 


G 126-2 


212907 


121059 


4500 


-1.42 


4.5 


0.9 


175 


32 


-14 


38 


-242 


1 


Y 


Sandagc86 


LP 873-63 


213139 


-214812 


5100 


-1 


4.5 


0.3 


55 


7 


-171 


79 


-119 


34 


Y 


Samir02 


WT 2218 


220744 


-221453 


6200 


-1.46 


4.5 


1 


-336 


232 


-275 


228 


-116 


162 


N 


Samir02 


LP 759-73 


221840 


-100812 


5400 


-1.85 


4.5 


0.9 


-277 


62 


-212 


29 


55 


51 


N 


Samir02 


LP 875-62 


222601 


-205031 


4600 


-1.03 


4.5 


0.4 


160 


57 


-90 


35 


20 


8 


Y 


Samir02 


HIP 112486 


224704 


295108 


5000 


-1.5 


4.5 


0.5 


-2045 


6852 


-1104 


2808 


-1412 


5086 


Y 


Samir02 


NLTT 55528 


230013 


215631 


5800 


-2.61 


4.5 


1 


-548 


351 


-307 


77 


-38 


89 


N 


Samir02 


G 273-76 


233457 


-203854 


4800 


-1.9 


4.5 


0.2 


36 


25 


-190 


61 


-26 


15 


Y 


Samir02 


G 157-88 


233757 


-74542 


5000 


-0.15 


4.5 


0.6 


250 


140 


-251 


127 


-23 


46 


Y 


Samir02 


G 273-105 


234117 


-194748 


4800 


-1.31 


4.5 


0.5 


200 


84 


-106 


48 


5 


4 


Y 


Samir02 


G 266-8 


235614 


-262559 


5200 


-0.81 


4.5 


0.3 


-106 


28 


-171 


42 


-16 


10 


Y 


Samir02 


G 30-31 


235754 


83646 


4200 


-1.86 


4.5 


0.4 


15 


13 


-71 


23 


-104 


18 


Y 


Sandagc86 


G 129-57° 


235831 


203424 


4700 


-0.82 


4.5 


0.3 


-45 


10 


-25 


3 


14 


2 


Y 


Sandagc86 



a hhmmss 
b ddmmss 

c Data taken with R=60000 



References. — (CLLA94) = Carney et al. (1994); (Hipp) = Hipparcos catalog; (NLTT) = Salim & Gould (2002); (Sandage86) = Sandage & Kowal (1986) 



Table 2. Comparison with literature 



This study Literature 



Star 


T of f 


[Fc/H] 


T ff 


[Fe 


/H] 


Source 


LP 644-47 


5400 


-1.47 




-2 


.07 


RN91 




5400 


-1.47 


5341 


-1 


.63 


CLLA94 


LP 705-89 


5800 


-0.97 




-1 


.12 


RN91 


G 72-18 


4900 


-1.74 


4982 


-1 


.86 


CLLA94 


G 39-1 


5000 


-0.06 


4867 


-0 


.61 


CLLA94 


HIP 19814 


5400 


-0.71 


5378 


-0 


.69 


CLLA94 


G 83-26 


4500 


-1.57 




-2 


.30 


RN91 


HIP 24289 


5000 


-0.66 


5032 


-0 


.77 


CLLA94 


G 249-37 


5000 


-0.13 


4954 


-0 


.27 


CLLA94 


LP 661-1 


5200 


-1.98 




-2 


.03 


RN91 


G 40-5 


5200 


0.06 


5012 


-0 


.07 


CLLA94 


G 9-13 


5000 


-0.58 


5049 


-0 


.27 


CLLA94 


G 43-30 


5300 


0.25 


5070 


-0 


.21 


CLLA94 


G 181-19 


4800 


-1.66 


4892 


-1 


.56 


CLLA94 


LP 813-13 


6000 


-1.42 




-2 


.10 


RN91 


LP 873-63 


5100 


-1.00 




-0 


.83 


RN91 


LP 759-73 


5400 


-1.85 




-1 


.63 


RN91 


G 273-105 


4800 


-1.31 






.17 


RN91 



References. — (CLLA94) = Carney ct al. (1994); (RN91) = 
Ryan & Norris (1991); 



This figure "figl.jpeg" is available in "jpeg" format from: 
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This figure "fig8.jpeg" is available in "jpeg" format from: 



http://arXiv.org/ps/astro-ph/0306076vl 



